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Spin exchange interaction between atoms in a spin-1 Bose-Einstein condensate causes atomic spin 
evolving periodically under the single spatial mode approximation in the mean field theory. By 
applying fast magnetic pulses according to a two-step or a four-step control protocol, we find ana¬ 
lytically that the spin dynamics is significantly suppressed for an arbitrary initial state. Numerical 
calculations under single mode approximation are carried out to confirm the validity and robustness 
of these protocols. This localization method can be readily utilized to improve the sensitivity of a 
magnetometer based on spin-1 Bose-Einstein condensates. 
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I. INTRODUCTION 

Spin exchange interaction between atoms in a spin- 
1 Bose-Einstein condensate (BEC) causes complex spin 
mixing dynamics and spin diffusion, which is a major 
obstacle to realize experimentally a high precision mag¬ 
netometer based on spinor BEC [l]-ll[- In order to im¬ 
prove the sensitivity of the magnetometer, a smaller spin 
exchange interaction is required, which may be imple¬ 
mented effectively by dynamical decoupling method us¬ 
ing optical Feshbach resonance techniques tana. In ad¬ 
dition, the small spin exchange interaction can be utilized 
to resolve the ambiguity of the spin texture in ferromag- 
netically interacting 87 Rb spin-1 BEC, where the spatial 
texture structure may be induced by the spin exchange 
interaction, the magnetic dipolar interaction, or both of 
them olsomi- 

However, a more experimentalist-friendly proposal to 
suppress the spin exchange interaction is employing the 
magnetic pulses and the microwave pulses, which are 
much easier to implement and tunc experimentally [Til 
EMI- By applying a magnetic field to an atomic spin-1 
BEC, only considered is the quadratic Zeeman effect <5 
which is proportional to the square of the field, because 
the linear Zeeman effect can be eliminated mathemat¬ 
ically by adopting the rotating reference frame, due to 
the conservation of the total magnetization of the spin- 
1 condensate [1 Ill EH- Under current experimental 
conditions, the effective quadratic Zeeman energy of ei¬ 
ther the magnetic field or the microwave driving field can 
be adjusted from -240 Hz to +240 Hz, which is about 
10 times larger than the spin exchange interaction for 
typical densities of a 87 Rb spin-1 condensate, ~ 10 14 
cm -3 E3, m, El . 

In this paper, we propose to localize the spin dynam¬ 
ics of a spin-1 BEC by periodically applying magnetic 
and/or microwave field pulses, which effectively suppress 
the spin exchange interaction. By applying two-step 
pulse cycles with positive 6 only, the condensate dynam¬ 
ics is localized if the relative phase of the initial state 
is close to zero; by applying four-step pulse cycles with 


both positive and negative 6 , the condensate dynamics 
is localized for an arbitrary initial state. The exploration 
of the robustness of the protocols shows that a wide pa¬ 
rameter regime exists for a spin-1 condensate under cur¬ 
rent experimental conditions. This proposal may find 
its potential application to improve the sensitivity of a 
practical high-resolution magnetometer based on spin-1 
BEC. 

The paper is organized as follows. In Sec. El we review 
the theoretical description of the free spin mixing dynam¬ 
ics under the single spatial mode approximation (SMA) 
in a spin-1 BEC in a magnetic field, whose quadratic 
Zeeman splitting 5 ranges from large negative values to 
large positive values. In Sec. mu we analytically design 
and numerically confirm the control protocols of mag¬ 
netic/microwave pulses to localize the condensate spin 
dynamics, where either a two-step or a four-step pulse 
cycle is employed. Furthermore, the robustness of the 
control protocols is explored in Sec. [IV] by assuming 5% 
random error of the pulse amplitude 5{t). Finally, a brief 
summary is presented in Sec. El 

II. FREE SPIN DYNAMICS IN A MAGNETIC 
FIELD 

Within the mean field theory, the free spin mixing dy¬ 
namics in a spin-1 BEC with either ferromagnetic or anti¬ 
ferromagnetic spin exchange interaction under the SMA 
in a magnetic field is described by the following equation 
of motion f|, 9s, 03 

Po = j-Po Vi 1 ~ Po) 2 ~ m 2 sin 9, 

6 = -^ + j( 1-2 Pq ) (1) 

ft a /(1 - po) 2 ^ m 2 

where c = C 2 -/V f dr\<j>{r )| 4 with N being the total number 
of atoms in the condensate and <f>{r) a normalized spatial 
mode function under the SMA, which is determined by a 
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scalar Gross-Pitaevskii equation with a spin-independent 
interaction, [— (ft * 2 * * * * * 8 9 /2M)V 2 + T4 xt (r)+Co|<)>| 2 ]<)>(r) = p4>(r) 
where M is the atomic mass and V is the external har¬ 
monic trapping potential. The spin independent coeffi¬ 
cient Co and spin exchange coefficient C 2 are given, respec¬ 
tively, by co = 47 rS 2 (ao + 202 ) /3M and C 2 = 47 iA 2 (a 2 — 
ao)/3M with the s-wave scattering length a 0 ( 0 - 2 ) for two 
spin -1 atoms in the compound symmetric channel of to¬ 
tal spin 0(2). For two popular ultracold spin-1 atomic 
gases in experiments, 87 Rb and 23 Na, c 0 |c 2 | is always 
satisfied and thus guarantees the validity of the SMA in 
most experimental situations [22| - [24| . 

The fractional population of spin component p a (ot = 
—1,0,+1) satisfies J2 a Pa = 1- The magnetization 
m = p+ — p- is a constant during the evolution, due 
to the isotropic nature of the spin exchange interac¬ 
tion. The relative phase among the three components 
is 9 = 6 + + 6 L — 20o with 9 a being the phase of the 
spin wave function. The quadratic Zeeman energy is 
8 = (E+ + E_ — 2Eq)/2 with E a the Zeeman energy 
shift of the component. In general, 8 « 72 B 2 Hz/G 2 
for 87 Rb BECs and 8 « 278 B 2 Hz/G 2 for 23 Na BECs, 
where the magnetic field B is in unit of Gauss. Due to the 
conservation of the magnetization to, the linear Zeeman 
energy (E- — E + )/2 can be eliminated mathematically 
by adopting a rotating reference frame. 

The total spin energy is a constant during the free 
evolution of the spin -1 condensate in a magnetic field 


e = cp 0 


(1 - po) + \/(l — Po ) 2 


to 2 cos(0) 


+ 5(1 - po). 


Starting from a given initial state, which is usually a 
ground state in a magnetic field in experiments, the con¬ 
densate evolves according to an iso-energy trajectory in 
the plane of po~9, by changing abruptly the magnetic field 
to a different value. By taking into account of the energy 
conservation, the Eq. m is further simplified as (p'o) 2 = 
(4/h 2 ){[£-8{l-po)}[(2cpo+8){l~po)-£\-{cp 0 m) 2 } thus 
the time evolution of po can be analytically expressed in 
terms of the Jacobian elliptic function cn(.,.) if S ^ 0 
and sinusoidal function if 8 = 0 [2i, jsj, (25| 


Po(t) = — [x 2 +X\- (x 2 - Xi) sin(7o + 2t\/2ce + c 2 to 2 )] 

( 2 ) 

for 8 = 0; 


p 0 (t) = x 2 + (x 3 - x 2 )cn 2 (70 + tyj2c8{xo - xi), k) 


(3) 


for c8 > 0; 

poit) = X2 - (x2 - xi)cn 2 ( 7 o + t^/—2c8(x3 - xi), k) 

(4) 

for c8 < 0. We have set h = 1. Here x\ < X 2 < 
X 3 (xi < X 2 for 8 = 0) are the roots of po = 0 
, k = ^ (X 3 — X 2 )/ (x 3 — xi) if c8 > 0, and k = 



0 / 7t 



FIG. 1. (Color online) (a) Typical trajectories in po-9 plane 
for <5/|c| = -0.5 (dashed line), -0.39,-0.1,0.2,0.7,1.5,1.59 
(solid lines), 1.8 (dotted line), from bottom to top. All the 
trajectories with solid lines evolve in a clockwise direction. 
The initial state (asterisk) is pot = 0.8 and 8i = 0. (b) De¬ 
pendence of the oscillation amplitude A (blue solid line) and 
the period T (green dashed line) of po on 8. The running phase 
modes corresponds to the region I and IV and the oscillatory 
modes corresponds to region II and III. The red dotted lines 
marked by Jo and J± denotes, respectively, the ground state 
quadratic Zeeman energy (the initial state coincides with the 
ground state and A is zero) and the resonant quadratic Zee- 
man energy (T is infinite and po —> 0 or 1). 


1 J (X 2 — xi)/(x 3 — xi) if c8 < 0 . 70 is determined by 

the initial state. Hereafter we assume c = — 1 thus the 

energy unit is |c| and the time unit is |c | -1 and to = 0 . 

Typical trajectories are illustrated in Fig. [3(a) for dif¬ 

ferent 8. Although starting from the same initial state, 

the trajectories could cover the whole po~9 plane by con¬ 

tinuously varying the quadratic Zeeman energy 8 from 
negative infinity to positive infinity. All the trajectories 
are classified into two modes: the oscillatory mode where 

8 is between [— w, 7r] and the running phase mode where 

9 goes beyond [— 7 r, 7 r]. As shown obviously in Fig. [3a), 

the oscillatory mode trajectories are evolving in a clock¬ 
wise (counterclockwise) direction if c < 0 (c > 0 ), while 
the running phase mode trajectories for large |5| |c| 

may take one of two opposite directions, depending on 
the sign of 8. This is a key point in order to localize the 
condensate spin dynamics. The boundaries between the 
oscillatory modes and the running phase modes satisfy 
one of the two requirements, poit) = 0 (8 = 5_) or 1 
(5 = 5+) if time is long. The corresponding period T be¬ 
comes infinite [see also Fig. [3b)] . Another special point 
8 = 8q denotes the coincidence of the initial state with 
the ground state thus the oscillation amplitude A is zero 
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FIG. 2. (Color online) (a) Typical trajectories under the 
modulation of 5(t). Each modulation cycle includes a free 
evolution with d = 0 (dashed line) and a controlled evolution 
with d/|c| = 1.5, 2.0, 3.0 (solid lines, from bottom to top). 
The initial state (red asterisk) is poi = 0.8 and 9i = — 0 . 17 T. 
(b) Dependence of the amplitude (blue solid line) and period 
(green dashed line) of po on d, the nonzero quadratic Zeeman 
splitting. Circles and diamonds are the period and ampli¬ 
tude calculated analytically with Eqs. 0 and m for large 
ds, respectively. 


but the period T is finite. 

The oscillation amplitudes and the periods are shown 
in Fig. mb). There are clearly four regions: (I) run¬ 
ning phase mode with increasing 9(t)', (II) oscillatory 
mode with 0 < po(t) < pod (HI) oscillatory mode with 
Poi < Po(t) < 1; (IV) running phase mode with decreas¬ 
ing 9(t). The amplitude of the oscillations A monoton- 
ically increases in regions (I) and (III) but decreases in 
regions (II) and (IV) with 5 increasing. The period of 
the oscillations T shows two resonant peaks at 5 = 5± 
where po(t) = 0 or 1 at long enough time. The period is 
almost a constant between these two peaks but decreases 
rapidly outside the peaks. Similar oscillation behaviors 
were observed also in antiferromagnetically interacting 
23 Na spin-1 condensates (c > 0) jl8, Us). 

We observe from Fig. 0a) that in the oscillatory mode 
9 increase or decrease with time if 9 ss 0 and po is around 
its extremes. We may utilize this property to localize the 
condensate dynamics around 9 ss 0 by canceling 9 in 
a period with 9 increasing (decreasing) during the first 
(second) part. For an arbitrary state, however, we may 
utilize both the oscillatory and the running phase modes 
to localize the dynamics since 9 may increase or decrease 
for a given po, depending on the value of <5. 


III. LOCALIZED SPIN DYNAMICS 


We consider first that the control period consists of two 
steps, a free evolution (d = 0) for a time slot t \ followed 
by an evolution in a magnetic field S = d for a time slot 
r 2 ■ We refer hereafter this protocol as two-step control. 
For a given initial state with 0i rs 0, it is easy to prove 
analytically that t\ depends on the initial state and t 2 
depends uniquely on d, which indicates that there is only 
one free parameter in the two-step control protocol. The 
time dependence of the magnetic field for the two-step 
control is 

= I °’ J ( T1 + T2 ) ^ t < j( Tl + T2 ) + Tl 

[ \ d, j{r i +t 2 ) + ti <t < ( j + l)(n + r 2 ) 


where j = 0,1,2, • • • is an integer denoting the number 
of control cycles. 

Typical controlled trajectories are illustrated in 
Fig. 0(a) for three values of d, where the initial state 
is p 0i = 0.8 and d; = —0.17T. We see clearly that the 
oscillations of both p 0 and 9 under the two-step control 
are smaller than that during free evolution, indicating 
that the condensate dynamics is indeed localized by the 
two-step control protocol. Starting from the same initial 
state, the lager the d is, the smaller the oscillation of po 
is. 

The condensate spin average is ( F ) = (F x )x + (F y )y + 
( F z )z for a state with (U 

(F x ) + i(F y ) = 2^/po(l - p 0 ) cos(d/2), 

(F z ) = 0. 

Once we localize po(t) and d(f), the condensate spin (F) 
is obviously localized. For a nonzero ( F z ), the localiza¬ 
tion occurs similarly. 

The cycle period T depends on the free evolution time 
Ti and the controlled evolution time r 2 . The free evo¬ 
lution time is determined by the evolution time of the 
system from its initial state poi and 9i to the symmetric 
state Pq{ti) = poi and 8(ti) = — 9i. In this way, the time 
ti is calculated analytically by using Eq. © 

Ti = (tt/2 — 7o)/\/2 ce + c 2 m 2 (5) 

In the limit of small 9 -C 1, ti « |di/[2c(l — 2poi)]| where 
we have used po(t) ~ p^. Similarly, the controlled evo¬ 
lution time t 2 is the evolution time of the system in the 
magnetic field 5 = d and can be calculated, by using the 
conditions po{T) = p$i and 9(T) = 9i, 

t 2 = \/ 27 o/\/-cd(a ;3 - x\) (6) 

In the limit of small 9 <C 1 and large d, t 2 « |dj/[2c(l — 
2poi) — d]|. In total, the cycle period is approximated as 


2c(l - 2p 0i ) [2c(l - 2p 0i ) - d] 

for small 9i and large d. We notice that T k, t\ for large 
d, as shown also in Fig. ©b). 
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We define the control oscillation amplitude as A = 
max(po)— min(jOo), which depends obviously on the initial 
state and the magnetic field d. The amplitude can be 
calculated analytically but is too lengthy to present here. 
In the limit of large d and small 0 *, the amplitude is 
approximately 


A — Ai + A2 


( 8 ) 


where 


A\ 

A 2 


POi(l POi ) „2 

4 ( 2 p 0 i — 1 ) *’ 

P 0 i(l — POi ) 

\2d/c\ - 4(2p 0 i - 1) 


We see that A approaches to A\ as d goes to infinity. 

In Fig[ 2 |b), we present the dependence of A and T on 
the control magnetic Held d. We see clearly that A and 
T decrease monotonically as d increases, manifesting the 
fact that better localization of the condensate dynamics 
is achieved in a higher magnetic field. We note that A 
and T approach their nonzero asymptotic values at large 
values of d. Actually, to reduce the oscillation amplitude 
A further down to zero, we have to employ the following 
four-step control protocol. 

We consider next that the control period consists of 
four steps, (i) a free evolution for a time ri, (ii) an evo¬ 
lution in a magnetic field with S = d\ for a time 72, (iii) 
a second free evolution for the time 73, and (iv) a sec¬ 
ond controlled evolution in another magnetic field with 
S = g? 2 for a time T4, as shown in Fig. EH a). We refer this 
protocol as four-step control. For simplicity but without 
loss of generality, we limit to the symmetric situations 
where d\ = —d 2 = d and T3 = Ti = r. It will be analyt¬ 
ically proved that T2 and 74 are uniquely determined by 
d and r. Thus there are only two free parameters, d and 
r, in the four-step control we considered here. 

It is straightforward to find the analytical solution to 
72 and 74, by using the initial state and the Eqs. m, 


t _ v^qo 

yj-cd{x 3 - Xl) 

v^qo 

T 4 — , = 

V cd(x ' 3 - x[) 

where K(k) is the elliptic integral of the first kind, q(, 
and x ' 3 4 are determined by po(r) and — 9 (r). We note 
here that the initial state for 72 is po(r) and 0 (t), and 
that for 74 is poi and — 6 i- The total period for a complete 
cycle is 


T = 2t 


\/ 2 qo 


+ 


y/H 


\J—cd{x 3 - Xl) -Jcd(x ' 3 - x[) 


( 9 ) 


Typical controlled evolution of the condensate is illus¬ 
trated in Fig. EJb), where the parameters are given in 
the caption. Compared to the two-step control protocol, 
there are two advantages. The first is that the initial 



FIG. 3. (Color online) (a) Schematic of a magnetic pulse se¬ 
quence. (b) Typical controlled trajectory of po(t) and 9(t) un¬ 
der a four-step pulse sequence of 5(t) for a four-step protocol. 
The red asterisk marks the initial state. The parameters are 
poi = 0.8, 0i = —0.27T, d/|c| = 10, and n = 7-3 = r = 0.0577. 
The values of 72 and 74 are calculated, 72 « 0.007277 and 
74 « 0.022677. (c) Amplitude of po under four-step pulse se¬ 
quences. Better localization of po (smaller A) is achieved for 
larger d and smaller t. 


state is arbitrary, particularly, 0 i goes beyond the small¬ 
ness requirement. The second is that the oscillation am¬ 
plitude and period approach to zero if d is large enough 
and t is short enough, as shown in Fig. He) and Eq. ®. 


IV. ROBUSTNESS OF THE CONTROL 
PROTOCOLS 

We have assumed the magnetic control pulses are per¬ 
fect in previous sections, but there are always uncontrol¬ 
lable errors in practical experiments, e.g., the microwave 
field 5 may have 5 % relative uncertainty [llj. Since the 
timing is pretty accurate in current experiments, we next 
evaluate the robustness of the two-step and four-step pro¬ 
tocols only under the 5 % uncertainty of 6 for many con¬ 
trol cycles. 

We define the fidelity of a protocol after many control 
cycles as 


F=\m f )\ 2 


( 10 ) 
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FIG. 4. (Color online) (a) Fidelity after 1 cycle (blue solid 
lines), 10 cycles (red dashed lines), 100 cycles (black dash- 
dotted lines), and 200 cycles (purple dotted lines) under the 
two-step protocol. For clarity, each curve is shifted up by 
0.01 from bottom to top. The insert shows the zoom-in view 
near the dip in the main panel, (b) Fidelity after 1 cycle (red 
lines) and 10 cycles (blue lines) under the four-step protocol 
for r/7r = 0.01 (dash-dotted lines), 0.03 (dashed lines), and 
0.05 (solid lines). The results show that both the two-step 
protocol and the four-step protocol are robust, i.e., the fidelity 
F is close to 1, in the presence of relatively 5% magnetic pulse 
errors. 

where |£i,/) is the initial and the final state of the spin-1 
condensate and satisfies |(£|£ )| 2 = 1. The state has three 
components, ||*) = (£+, £o, £-) T with 
and p a and 9 a being the fraction and the phase of the 
component a, respectively. The fidelity measures how 
close the initial and the final states are. The fidelity is 1 
for ideal pulses but lower than 1 in the presence of pulse 
errors. The larger the errors are, the lower the fidelity is. 
Higher fidelity indicates more robustness of the protocol 
to pulse errors. 

We assume the magnetic field error is distributed with 


equal probability in the range [0.95,1.05]d with an aver¬ 
age of d. For the four-step protocols, the errors for d\ 
and g ?2 are independent. We numerically calculate the 
dependence of the fidelity F on d and show the results 
in Fig.[4ja) for two-step protocols and[4](b) for four-step 
protocols. 

As shown in Fig. Oja), the fidelity is above 99% for 
most d , except a special dip near d/\c\ « 1.2 where the 
period T is most sensitive to the change of d [i.e., the 
largest derivative of T with respect to d in Fig. [2jb)] . 
This result manifests that the two-step control protocols 
are pretty robust under the uncertainty of the magnetic 
field, if we choose a field away from the dip. 

For the four-step protocols, as shown in Fig. |U[b), we 
find that the fidelity is very close to 1, though it decreases 
as d increasing or r decreasing. By taking into account of 
the requirement of small r and large d to better localize 
the condensate dynamics, a delicate balance between the 
localization and the robustness is required in practical 
experiments. 


V. CONCLUSION 

We propose to to localize the spin mixing dynamics in 
a spin-1 Bose condensate by periodically applying mag¬ 
netic pulse sequences, according to the two-step proto¬ 
col for an initial state with small initial relative phase 
or the four-step protocol for an arbitrary initial state. 
Numerical calculations confirm the validity of the pro¬ 
posal for a ferromagnetically interacting spin-1 conden¬ 
sate under the single spatial mode approximation. We 
further illustrate the robustness of the localization pro¬ 
tocol with numerical calculations by assuming 5% un¬ 
certainty of the magnetic pulse amplitude, which might 
occur in practical experiments [l8| . Our proposal may be 
utilized to realize hig her precision magnetometers based 
on spinor BEC [4|,[26j or to explore the weak dipolar inter¬ 
action effects in 87 Rb spin-1 condensates by suppressing 
the spin dynamics induced by the spin exchange interac¬ 
tion UMMMm. 
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